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TECHNICAL NOTE ^l66 


AN EXPERIMENTAL AND THEORETICAL STUDY OF THE EFFECT 
OF FUEL ON PITCHING-TRANSLATION FLUTTER 
By John L. Sewall 


SUMMARY 


Analytical flutter studies were made for 2 two-dimensional fuel- 
loaded wing models, and the results are compared with experimental results 
for hending-to- tor sion frequency ratios near 1. One of the models was 
made so that water, simulating fuel, could he carried intern ally in three 
compartments separated from each other hy sealed spanwise partitions. 

The flutter speeds of this model for all fuel loads were highest for the 
conrpartmental-enrptying sequence proceeding from front to rear. In the 
other model, fluid was carried externally in a geometrically scaled 
standard airplane fuel tank that was pylon mounted a distance of about 

2^ times the tank radius beneath the wing. Experimental results for this 
model have been reported in NACA Research Memorandum L55F10. 

The results of flutter-speed calculations agreed well with experi- 
mental flutter speeds and flutter-speed trends when the analysis employed 
effective values for the mass and the mass moment of Inertia together 
with theoretical slender -body air forces on the external tank and when 
the combined structural and fluid damping was considered to be zero. 

With the introduction of damping this agreement was improved for the 
Internal- tank configuration and made worse for the external- tank 
configuration . 


INTRODUCTION 


The flutter of airpla n e wings with heavy fuel loads inherently 
involves the superposition and interaction of two highly complex dynamic 
phenomena, one associated with wing flutter and the other with fuel 
motion. This dual nature of the problem has motivated two main avenues 
of inquiry, one of which has been concerned with the fuel motion itself 
and the other with the effect of fuel motion on flutter. Fuel dynamic 
studies, such as references 1 to 6, have been chiefly concerned with 
finding ways to represent or approximate the complicated motion of the 
fuel in oscillating tanks. 
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Investigations of the effect of fuel motion on flutter have dealt 
with either realistic configurations, as, for example, in reference 7, 
or with simple models, as is the case in reference 8 and in the present 
paper. The purpose of using simple models for an investigation of this 
kind is to isolate as nearly as possible the effects of fuel motion in 
order to obtain results applicable to realistic configurations but not 
appreciably influenced by certain structural and aerodynamic effects 
peculiar to those configurations. With this purpose in mind, a two- 
dimensional wing apparatus was used in reference 8 to obtain flutter 
data for an external wing-tank configuration in which the tank was pylon 
mounted beneath a semirigid untapered wing mounted on bending and tor- 
sional springs located outside the air stream. 

This same objective applies to the present paper which is primarily 
concerned with the correlation of flutter theory with experiment not 
only for some of the data presented in reference 8 but also for data 
obtained by using a two-dimensional internal wing- tank configuration. 

The two-dimensional features of the apparatus and models were particu- 
larly well suited to the simplest application of flutter theory because 
the assumption of two-dimensional air flow over the wing was more fully 
satisfied and because the structural behavior of the system could be 
adequately represented by just two vibration modes, namely, vertical 
translation (simulating first bending) and pitch (simulating first 
torsion). The mass and mass moment of inertia of the fuel were approxi- 
mated by means of effective values which were determined during simple 
harmonic oscillations of the wing-tank system in vertical translation 
and pitch, respectively. As shown in references 4 to 6 , these values 
are usually less than the mass and mass moment of inertia of the same 
volume of fuel at rest (or considered frozen). The results of refer- 
ence 8 demonstrated the usefulness of references 4 to 6 by showing that, 
when the fluid in the tank was replaced by Bolid bodies having masses 
and mass moments of inertia equal to the corresponding effective values 
for the fluid, the model fluttered at the same speeds as it did with 
fluid. The present paper seeks to determine how well flutter speeds can 
be predicted by classical flexure-torsion theory based both on effective 
inertial properties of the fluid and on damping of the fluid according 
to a commonly used approximation. 

The flutter analysis used is the same as that given in reference 9 
except that the mass In translation is different from the mass in pitch. 

In addition, the aerodynamic forces on the external tank are included. 

The effect of these additional forces on flutter and divergence is 
explored for the external-tank model as a secondary objective of this 
paper. An appendix contains the flutter and divergence equations together 
with expressions for the tank aerodynamic coefficients derived by appli- 
cation of s lender-body theory. 
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SYMBOLS 

dimensionless distance of elastic axis relative to midchord, 
positive rearward 

rearward fuel compartment of model A 

semichord of wing, ft 

center fuel compartment of model A 

distance from wing quarter-chord to elastic axis, positive 
rearward, ft 

distance from wing quarter-chord to volumetric centroid of 
external tank, -positive rearward, ft 

forward fuel compartment of model A unless otherwise noted 

frequency, cps 

uncoupled translation frequency, cps 
uncoupled pitching frequency, cps 

value of damping coefficient which BatisfieB flutter deter- 
minant as a fu nc tion of V/tx% in flutter analysis (see 
fig. 12) 

experimentally determined damping coefficients in translation 
and pitch, respectively 

vertical translation of elastic axis of wing section, or 
bending, positive down 

mass moment of inertia of model in torsion about elastic 
axis, ft-lb-sec^ 

effective mass moment of Inertia of fluid about elastic axis, 
ft-lb-sec^ 

mass moment of inertia of fluid considered as a solid about 

O 

elastic axis, ft-lb-sec 
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volume moment of inertia of external t ank of model B about 
elastic axis, ft^ 

reduced frequency, bcc/V 

spring constant of model in translation, lb/ft 
spring constant of model in pitch, ft-lb/radian 

outside radius of external t ank of model B, ft 
dimensionless radius of gyration of model about elastic axis 

i ^ 2 

wing semispan, ft 
airstream velocity, fps 

experimental flutter speed of model in empty condition, used 
as a reference speed in figs. 6(b), 7(b), 8(b), 9, and 11, 
fps 

volume of external tank of model B, ft3 
mass of model in torsion, lb-sec2/ft 
mass of model in bending, lb-sec 2 /ft 
effective mass of fluid, lb-sec 2 /ft 

actual mass of fluid (considered aB a solid), lb-sec 2 /ft 

location of elastic axis of model from leading edge, percent 
chord (see table I) 

location of model center of gravity from leading edge, per- 
cent chord (see table I) 

location of model center of gravity relative to semichord, 
measured from elastic axis, positive rearward 

angle of pitch of wing section about elastic axis, positive 
leading edge up 
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k ratio of mass of cylinder of air of radius b to mass of 

model, Sipb^s/m 

0 ) angular frequency, rad Ians /sec 

p air density, Ib-sec^/ft^ 

Subscripts : 

d torsional divergence 

f flutter 

h translational degree of freedom 

0 tank-empty condition 

T external tank of model B 

a pitching degree of freedom. 

When used together the letters R, C, and F denote a particular 
fuel-emptying sequence for model Aj for example, RCF identifies the 
sequence in which the rearward compartment is assumed to be emptied first, 
then the center compartment, and finally the forward compartment . 

EXPERIMENTAL INVESTIGATION 

Models and Apparatus 

The experimental data presented in this paper were obtained by the 
use of two simplified wing models, cross-sectional views of which are 
shown in figure 1. One of these models, designated model A in figure l(a) 
was made so that fluid could be carried internally in three compartments 
which were separated from each other by sealed spanwise partitions. In 
the other model, designated model B, the fluid was carried externally in 
a geometrically scaled standard airplane fuel tank that was pylon mounted 

a distance of about 2^ times the tank radius beneath the wing as shown in 

figure l(b) . Both wings were of aluminum-alloy construction except far a 
wooden trailing edge on model A and an Internal hardwood core formed to . 
the airfoil contour of model B. The skin thickness for each model was 
0.032 inch, and the thickness noted for model A In figure 1(a) includes 
a layer of fiber gla ss and Paraplex which was wrapped around the outside 
of the wing. The center of the tank pylon of model B was located above 
the center of the tank midsection and was fastened to the wing through, a 
bracket which permitted chordwise variation of the center -of -gravity 
position of the pylon-tank assembly. Tables I and H list some 
additional data for the models. 
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Both models were fluttered In the same two-dimensional apparatus 
described in reference 8. A schematic view of this apparatus is shewn 
in figure 2. Arrows indicate the translational and pit ching degrees of 
free dan restrained by bending (translational) and torsional (pitching) 
springs which were mounted outside the test section of the Langley 2- 
by 4-foot flutter research tunnel. 

The spring constants of one set of bending and torsional springs 
were selected so that the uncoupled-wing bending-to-torsion frequency 
ratios for both fuel-loaded models approached and exceeded 1 with increasing 
tank fullness. Values for the spring constants are given in table I. A 
second set of spring constants listed in the table for model B yielded fre- 
quency ratios greater than 1 for flutter .conditions of this model without 
fuel over a range of forward and rearward center-of -gravity positions of 
the pylon-tank assembly from 30 to 60 percent wing chord. 

Both models were fluttered in air at atmospheric pressure with fuel 
loads simulated by water and by solid weights. In model A the solid 
weights were mounted on the pitch axis outside the test section and in 
model B they were mounted inside the tank. 


Effective Inertial Properties of Fuel-Loaded Wings 

The apparatus was equipped with restraints so that the model could 
be excited in one degree of freedom independent of the other. This was 
an important feature because it permitted the measurement of uncoupled 
frequencies in either pure translation or pure pitch. These frequencies 
were used to determine the amount of fluid effective as mass and mass 
moment of inertia during flutter by application of essentially the same 
technique as that used in references 4 to 6. In order to determine the 
best approximations to effective values of these inertial properties at 
flutter, the frequencies read on the vibration records covered the same 
range of amplitudes encountered during flutter. 

It is Intended by repeated references to- effective values to empha- 
size the fact that the effective mass and mass moment of inertia of an 
oscillating tank containing fluid are different from the mass and mass 
moment of inertia of the tank contai ning the same volume of fluid at rest, 
or in a solid condition. This behavior, of course, occurs because the 
fluid, during the oscillation, does not follow the motion of the tank 
perfectly. Indeed, as previously noted, references 4 to 6 have shown 
effective values of these properties to be generally less than the cor- 
responding solid values, as might be intuitively expected. 

The center-of -gravity variations with tank fullness are shown in 
figure 3. These variations are based on experimental values obtained 
with fluid in the tanks for the various fullnesses and fluid distributions 
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• represented by the symbols. The straight lines connecting these symbols 
indicate that 3 1 nM.r interpolation was employed to obtain centers of 
gravity for intermediate fullnesses and fluid distributions. For model A 

. the symbol explanation in figure 3(a) applies throughout the paper and 
represents the fuel distribution in increments of full compartments for 
the six fuel-emptying sequences simulated in the experimental flutter 
program. In subsequent figures pertaining to model A the symbol explana- 
tion is modified to represent intermediate half -full compartments for 
which flutter tests were also conducted in the experimental, program. 

Each sequence is represented by the order of the letters F, C, and R 
to indicate the order in which the compartments would be emptied during 
flight. For example, FCR applies to the sequence in which the forward 
compartment would be emptied first, then the center compartment, and 
finally the rearward compartment. 

Figure 4 shows the variation of effective mass with tank fullness 
and figure 5* the variation of effective mass moment of inertia with tank 
fullness. The data in both figures have been referred to the mass and 
mass moment of inertia of the tank-empty condition for both models. 
Numerical values beside the points on the dotted curves and lines give 
the amount of fluid effective as mass and mass moment of inertia during 

• translatory and pitching oscillations, respectively. In obtaining the 
data in figure 4, the variation of the measured translational frequency 
with mass was first established for several different solid weights 
mounted on each model in the empty condition. The effective mass of the 
model containing a given amount of fluid is, then, the empty mass plus 
the equivalent solid mass corresponding to the translational frequency 
measured for that fluid load. For model A the effective mass and actual 
mass were Identical in all except the two fuel loadings Indicated in 
figure 4(a) . 

The effective mass moments of inertia presented in figure 5(a) and 
table III were determined from the relation 


i 

(a*a f 

Values of I a g greater than 1 cannot be readily explained but may 

be due to the swirling action of the fluid in a partly filled compartment. 
Moreover, where effective mass moments of inertia are only slightly greater 
than the mass moment of inertia of the empty model ^that is, l^jl^ Q ** 1^ , 

determination of the fluid-inertia ratio :Ls su ^ J ec t to small- 


difference errors 
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The method of determining the effective mass moments of Inertia pre- 
sented In figure 5(b) is similar to that used for the data in figure 4. 
Solid weights, mounted inside the empty external tank, were designed to 
the solid mass moments of inertia of the fuel for tank fullnesses up to 
100 percent, and the variation of measured pitching frequency with these 
solid inertias was established. The effective mass moment of inertia of 
the model containing a given amount of fluid was then found by adding the 
model-empty mass moment of inertia to the equivalent solid mass moment of 
inertia corresponding to the measured pitching frequency for that fluid 
load. 

Experimental Results 

Model A .- As previously mentioned, the internal-tank configuration 
was tested with the wing filled with water in increments of half -full 
compartments for all six of the fuel-emptying sequences. Vibration and 
flutter data obtained in this program are given in table IV(a) . The 
damping coefficients are based on logarithmic decrements of the amplitude 
decays for the natural frequencies listed in the table. 

The damping coefficients and flutter speeds listed in table IV(a) 
are shown for the complete range of tank f ulln ess in figures 6 to 8 , 
in which the flutter speed is referred to the tank-empty flutter 
speed V f Q . Figures 6 (b), 7(b), and 8 (b), together with figure 3(a), 

show that, except for the emptying of the center compartment in the RFC 
and FRC sequences (fig. 8 (b)), the flutter speed consistently increased 
as the center of gravity was shifted rearward. If no attention is given 
to translation-to-pitch frequency ratios (equivalent to bending-to- 
torsion ratios), this trend appears at first glance to conflict with the 
established trend toward reduced flutter speed due to rearward center-of- 
gravity shifts for low frequency ratios. However, it has been shown in 
reference ID that, when the frequency ratio lies within the range of 0.8 

to about 1.3 for a = -0.4, r a 2 = 0 . 25 , and" 1 /k = 20 , the flutter speed 

experiences a marked increase as the center of gravity is moved back. 

This effect is s imilar to that found for the present configuration in 
which a = -0.4, r a ^ ranged from 0.17 to 0.3^5, l/* ranged from 31 

to 87, and for which the range of frequency ratio was from 0.76 to 1 . 06 . 

Further examination of the experimental flutter data in figures 6 (b), 
7 (b), and 8 (b) reveals the existence of an optimum fuel- emptying sequence, 
that is, the compartmental sequence resulting in the consistently highest 
flutter speeds for all fuel loads . This is. the sequence which gives the 
most rearward center-of-gravity shift for any tank fullness and is des- 
ignated FOR in figure 6 (b), which also shows that the flutter speed for 
this sequence was as much as 22 percent higher than it was for the opposite 
sequence (RCF). 
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Model B .- The experimental results reported in the present paper 
supplement those reported in reference 8 for the uhbaffled external-tank 
configuration. These results are listed in table 17(b). The damping 
coefficients for all partially full conditions were determined from 
frequency-response curves obtained at zero airspeed and correspond to 
the translational and pitching flutter amplitudes also listed in the 
table. The damping coefficients for the tank-empty and tank-full con- 
ditions were determined in the same manner as were those for model A. 

The flutter speeds listed in table 17(b) are shown, referred to the 
tank-empty flutter speed, in figure 9 where they are given as functions of 
the tank f ullne ss. 7alues of the frequency ratio are shown beside 

each of the data points. The significance of these values with regard to 
fluid damping is discussed in reference 8 and in a subsequent section of 
the present paper. This strong frequency-ratio effect on flutter speed is 
not only of practical importance in realistic configurations but has also 
been encountered in such configurations as, for example, in the case of 
the dynamically scaled wing-tip tank model of reference 7* 

The limiting effect of fluid turbulence on flutter amplitudes has 
also been discussed in reference 8 and is examined further in the present 
paper. As table 17(b) shows, the model fluttered over a wide Bpeed range 
and a correspondingly limited amplitude range for all partially full con- 
ditions. Notice, however, that the flutter amplitudes consistently 
Increased with further increase of speed for tank fullnesses ranging from 
25 to 75 percent; whereas, for greater fullnesses the amplitudes tended 
to reach a peak and thereafter Consistently decreased with increasing 
Bpeed. This behavior may be seen more clearly in figure 10 in which the 
translational and pitching components of the flutter amplitudes are shown 
as functions of airstream velocity for the 85-, 90 -, and 95-percent full 
conditions. 7alues of the damping coefficients g^ and g a corre- 
sponding to the flutter amplitudes are shown beside each point. As noted 
in table 17(b) with the tank 95 -percent full, the model essentially 
stopped fluttering when the pitching amplitude became 0° and the trans- 
lational amplitude was barely visible. Thus, there appears to be a closed 
or bounded flutter-speed region for this case. On the basiB of this 
experience and of the pitching-amplitude trends for the other two fuel 
loads, it was possible by extrapolation (indicated by the dashed curve) 
to establish closed flutter regions for these cases also. The upper 
flutter-speed boundary determined in this manner is shown in figure 9 
in terms of the tank-empty flutter speed and can be seen to be anywhere 
from 62 to 97 percent higher than the lower flutter-speed boundary. 

In addition to the fuel-load studies, flutter tests were also con- 
ducted for the empty external tank located at a few other chordwise center- 
of -gravity positions. A different set of bending and torsional .springs 
was vised, and the frequency ratios obtained were greater than 1. The 
results of these tests are given in table 17(c) and are shown as a function 
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of 3^ in figure 11 in which the flutter speeds are referred to the 
same tank-empty flutter speed used for figure 9* The range of values 

corresponds to a range of tank center-of -gravity positions from 30 percent 
wing chord to 60 percent wing chord. The dangling coefficients are based 
on logarithmic decrements of the amplitude decays corresponding to the 
natural frequencies listed in the table. As was expected, the range of 
forward center- of -gravity positions at which flutter could be obtained 
was limited by the divergence-speed boundary. 


ANALYTICAL INVESTIGATION 


Method and Application of Analysis 

The two-dimensional charact eristics of the models and apparatus, as 
illustrated in figure 2, made the simplest form of flexure -tors ion flutter 
theory, as given by the equations in the appendix, particularly applicable 
in the present paper. For one thing, the flutter mode could be adequately 
represented by Just two degrees of freedom, namely, rigid vertical trans- 
lation (equivalent to bending) and rigid pitch (equivalent to torBion) . 
Moreover, the oscillatory aerodynamic forces and moments acting on the 
wing could be more truly reproduced by use of the classic incompressible - 
flow coefficients of reference 9* 

The dynamic behavior of the fuel was represented in the flutter 
analysis by means of the effective values of mass and mass moment of 
inertia given in figures k and 5 and the measured damping coefficients 
given in table IV. The validity of using effective values for the inertial 
properties has been established experimentally in reference 8 wherein the 
flutter speeds of model B containing fluid were generally duplicated when 
the fluid was replaced by solid weights equivalent in effective mass and 
mass moment of inertia to the fluid. The measured damping coefficients 
were used to approximate the effect of damping at flutter by application 
of the well-known assumption given in reference 11, by which the damping 
coefficient g h or is considered to be indicative of the ability of a 

given configuration to change the flutter speed from that corresponding 
to g = 0, as may be seen, for example, in figure 12 which will be dis- 
cussed later. The merits of this assunption are examined separately for 
each model. 

Oscillating aerodynamic forces and moments acting on the external 
tank of model B were taken into account by application of the s lender- 
body theory to obtain expressions from which aerodynamic coefficients for 
the tank were derived. As shown In the appendix, these coefficients were 
then sinply added to the two-dimensional aerodynamic coefficients in the 
flutter analysis. 
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** Analytical Results 

# The results of the analytical program are given in table V in the 

same order as the experimental results which are given in table 17 . For 
model A, damping was introduced into the calculations according to the 
assumption specified in the preceding section, and the amount of damping 
assumed to be present was arbitrarily based on the larger value of damping 
coefficient listed for each fuel load in table IV(a) . This predominant 
damp ing coefficient , either g^ or g^, is repeated for each "case in 

table V(a) and is plotted as a function of tank fullness and fuel distri- 
bution in figures 6(a), 7 (a), and 8(a). The lower boundary of the shaded 
region in figures 6(b), 7(b), and 8(b) represents the calculated flutter 
speed with damping neglected, and the upper boundary represents the cal- 
culated flutter speed with damping included. As can be seen, the intro- 
duction of damping resulted in flutter-speed increases of from 1 
to lk percent. 

Hie analy tical results for model B are listed in tables V(b) 
and V(c). Since the tank was externally mounted, particular attention 
was directed to the effect of including oscillating aerodynamic forces 
and moments acting on the tank. The effect of damp in g is considered in 
the correlation of theory and experiment . For the fuel-loaded tank 
fixed at 40 percent wing chord, the linear superposition of two-dimensional 

« wing forces and slander-body tank forces, according to the preceding sec- 

tion, reduced the flutter speed by 7 to 11 percent. For the tank-eagjty 
configuration, this reduction was increased from 7 percent for a frequency 
ratio of O.67 to ID percent for a frequency ratio of 1.13. When the empty 
tank was moved rearward to 60 percent of the wing chord, the inclusion of 
tank forces resulted in a reversed effect; that is, the flutter speed was 
increased by 10 percent for very nearly the same frequency ratio. (See 
also fig. 11.) 

As a matter of interest, stimulated by the fact that model B without 
fuel encountered divergence rather than flutter with the tank center of 
gravity at 30 percent wing chord, the effect of tank aerodynamic forces 
on the divergence speed was examined by using the divergence -speed equa- 
tion given in the appendix. The results are presented in table V(c) and 
are also shown in figure 11 . As can be seen, the calculated divergence 
speeds obtained by including tank forces are approximately 10 percent 
lower than those obtained by neglecting tank forces. 


CORRELATION OF THEORY AND EXPERIMENT 


General Discussion 


The exper ime ntal and analytical results of this investigation are 
compared in figures 6 to 8 for model A and in figures 9 , 11 , and 12 for 
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model 5. As can be seen, the use of effective values of mass and mass 
moment of inertia in flutter calculations gave answers which agreed satis- 
factorily with experiment for nearly all fuel conditions. This result has 
also been obtained experimentally for model B in reference 8 which shows 
excellent agreement between the flutter speeds of the model with fluid and 
the flutter speeds of the model with solid weights equivalent in effective 
mass and effective mass moment of inertia to the fluid. For model A the 
calculated flutter speeds adequately follow the experimental trends for 
all fuel-emptying sequences and are consistently conservative, or less than 
the experimental flutter speeds, by margins ranging from 2 to 20 percent. 
When fluid damping was taken into account in the manner previously 
described, these margins narrowed to less than 10 percent, and the calcu- 
lated flutter speeds were still conservative. 

Neglect of aerodynamic forces and moments on the external tank of 
model B led to calculated flutter speeds that were as much as 11 percent 
unconservative, or greater than the experimental flutter speeds, for the 
fuel-loaded configuration. (See fig. 9 for ‘the 25 -percent-full condition.) 
With tank forces included, the calculations became generally conservative. 
As can be seen, theory appears to follow closely the experimental trend up 
to the 95 -percent -full condition for which the calculated flutter speed was 
infinite, whereas the experimental flutter speed was not. 

Calculated flutter speeds were also conservative for model B in the 
empty condition with the tank center of gravity located at various chord- 
wise positions ranging from 50 percent wing chord to 60 percent wing chord. 
As is shown in figure 11, the difference between calculated and experimen- 
tal flutter speed amounts to as much as 21 percent. Furthermore, theory 
including tank aerodynamic forces exhibits a somewhat more consistent 
tendency to parallel experiment than does theory neglecting tank forces. 

It can also be noted that inclusion of tank forces resulted in a calcu- 
lated divergence speed which agreed well with the experimental divergence 
speed for the tank at 50 percent wing chord but was 6 percent less than 
the experimental flutter speed when the tank was at 4-0 percent wing chord. 


Some Further Remarks on Damping. 

The Improved agreement between theory and experiment due to consid- 
ering g equal to the experimental value of g^ or g^ whichever is 

larger, has already been discussed for model A. For model B, however, 
the situation is different. Figure 12 shows g as~a function of flutter- 
speed coefficient V/bci^ for all the tank fullnesses tested with model B. 

The solid curves define theoretical borderline conditions between stable 
and unstable motion over a limited range of . g. In order to be consistent 
with the assumption regarding damping, the experimental flutter results 
from table IV(b) are positioned in figure 12 according to experimental 
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values of and g^. The points identified by the open circles repre- 
sent speeds at -which the borderline conditions between stable and unstable 
motion were first encountered experimentally. The solid circles represent 
experimental flutter speeds within the unstable region for all partially 
f ull conditions in which the flutter amplitudes were limited by fluid tur- 
bulence. The theoretical solutions at. g = 0 correspond to the calcu- 
lated flutter-speed ratios presented in figure 9 with tank aerodynamic 
forces and moments included. The relation between theory and experiment 
for the tank-empty configuration (fig. 12(a)) is similar to that found 
for all cases for model A. 

It is important to note that the effect of dampin g appears to be 
more pronounced at the higher fuel loads as is indicated by the consistent 
decrease in the slope of the theoretical curve . This behavior is attrib- 
uted pr imari ly to the fact that the frequency ratio approached and 
exceeded 1 as the tank fullness approached 100 percent. Hie dashed line 
shown in figure 12(g) represents the stable condition (i.e. , no flutter) 
shown in figure 10 for a torsional component of the flutter amplitude 
equal to zero. 

This sensitivity of theory to 1 amounts of damping was also 
typical of the cases of the empty-tank configuration tested with the tank 
in different chordwise positions. Inclusion of damping according to the 
assumption g = g^, where g^ is less than g^ resulted in flutter 

speeds that were unconservative by as much as 78 percent with the tank 
center of gravity at 40 percent wing chord. 


CONCLUSIONS 


This paper reports experimental and analytical investigations of 
flutter of two simplified fuel-loaded wing configurations, the two- 
dimensional characteristics of which are particularly suited to the 
simplest application of classical flexure-torsion flutter theory. The 
wings of both models were semirigid, untapered, and mounted on bending 
and torsional springs located outside the air stream. Water was used to 
simulate fuel carried internally in one configuration in three compart- 
ments separated from each other by sealed spanwise partitions and car- 
ried externally in the other configuration in a geometrically scaled, 
pylon-mounted, standard airplane fuel tank. The results of the investi- 
gations appear to justify the following conclusions for these 
configurations : 

1. For the wing with fuel stored in compartments within the wing 
structure and with bending-to-torsion frequency ratios ranging from O.76 
to 1.06, opt imum (i.e., consistently highest) flutter speeds were found 
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for the campartmental— emptying sequence proceeding from front to rear and 
corresponding to the most rearward shift in the center of gravity. It 
be noted that this trend is in contrast to the established trend 
toward reduced flutter speeds due to rearward center-of -gravity shifts 
for low frequency ratios. 

2. The conclusion of NACA Research Memorandum L55F10 is herein 

reaffirmed, namely, that effective or average values of mass and mass 
moment of inertia should be used in flutter calculations of fuel-loaded 
wings . ’ 

3. For the external wing-tank configuration the inclusion of s lender- 
body approx imat ions to the aerodynamic forces and moments acting on the 
external tank at flutter tended to improve the agreement between flutter 
theory and experiment. 

4. Arbitrarily assuming the damping coefficient at flutter to be 
either the coefficient in translation or pitch gave conservative answers 
that agreed well with experiment for the internal-tank configuration but, 
for the wing with the external fuel tank, flutter speeds calculated on 
the basis of either damping coefficient were unconservative by wide 
margins. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., August 20, 1957- 
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APPENDIX 


FLUTTER AND DIVERGENCE EQUATIONS FOR TWO-DIMENSIONAL 
WING-TANK CONFIGURATION 


Flatter 

The calculated flutter speeds presented in the body of this paper 
were obtained by use of the equations given in this appendix. As pre- 
viously noted, the aerodynamic forces and moments acting on the exter- 
nal wing-tank configuration of model B were approximated by a simple 
linear superposition of two-dimensional and slender-body aerodynamic 
theories. The two-dimensional, forces for the wing are those given In 
reference 9 , and the slender-body forces acting on the external tank are 
essentially the same as those given In reference 12 far closed bodies of 
revolution. The fineness ratio of the external tank of model B, namely 
7-15 (see table I), Is, as shown in reference 13, within the range of 
fineness ratio for -which the basic slender-body assumption of two- 
dimensional incompressible flow in planes normal to the free-stream 
direction is applicable. 

Linearly superposing wing and tank aerodynamic forces and moments 
by simply adding them together results In the following equations of 
motion: 


-m'cD 2 !! + k^ (l + igjjh - mba^p^a - (p + Pip) = 0 


( 1 ) 


— Sbx^ + k a (l + lgja - - (% + M a ^ ip ) = 0 


( 2 ) 


in which the mass in translation m* is greater than the mass in pitch m 
because It includes the torsional springs (see fig. 2) and where 


Prp = pvn> 2 h - pvu) 2 a^l a - lip + i 


) 


- *a + 1 |) + + ^ v ) 
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In these equations, x is the chordwise coordinate measured from the 
wing quarter-chord, positive rearward; y is the spanwlse coordinate; 
m 1 is the mass of the model in translation per unit span; m is the mass 
of the model in pitch per unit span; I a is the mass moment of inertia of 

the model in pitch about the elastic axis, per unit span; P is the aero- 
dynamic force acting on the wing per unit span; and is the aerody- 

namic moment acting on the wing about the elastic axis per unit span. The 
foregoing integrals for the tank were evaluated numerically over the tank 
length on the basis of the tank ordinates given in table II, and the 
results are given in table II (b) for all tank positions considered in 
the body of the paper. 

Introducing uncoupled frequencies by means of the expressions 
o^ 2 = It^/m' and o^ 2 = k^I a performing certain algebraic manipula- 

tions leads to the equilibrium equations 


t A n + ° 0 . 

(3) 

| + aE^j. = 0 

oo 


from which is obtained the determinants! flutter equation 
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A 11 B 11 


D n E n 




A ch + A ch, 


.5sL. (. 

K * V 


A ca, + A c<x,t) 


D li - T? - (\h + **h,x) 

- f 1 - (if) 2 ! 1 + *«■)] ^ ' ( A *a + Vt) 


2rtp"b 2 s 


2xpb^e 


2 _ 1 a 


x a ^ K £' A a mb 2 

aaie aerodynamic coefficients for the wing, frcm reference 14, are 


. n 2Q . ,2F 

A ch ■ _1 ” E" 1- h 


Aca = a + 


5 .(|. a ) f+ 1 [ i+ g + (l- a )f] 

Aah | “ (| + a ) A <* 


A aa = -^ 
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-where F and G are the well-known aerodynamic functions derived In 
reference 9 , and the aerodynamic coefficients for the tank are 



Divergence 

Since torsional divergence is a static instability associated with 
the torsional degree of freedom, the speed at which an external wing-tank 
configuration will diverge can be obtained by setting h = 0 (eq. (l)) 
and cd = 0, whence F + iG -*1. Under these conditions the divergence 
speed v d , which can be obtained frcsa a reduced form of equation (2) , is 
given by 


v d 

Mi 



( 6 ) 


The divergence speeds presented in the main body of this paper were cal- 
culated by means of this equation. 
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TABLE I.- PERTINENT MODEL CHARACTERISTICS 


Model A 


b, ft 

2s, ft 

Airfoil section NACA 

Fineness ratio of external tank 

Maximum fuel length (chardwise), in.: 

Forward compartment ............. 

Center compartment .... 

Rearward compartment 

Maximum fuel depth, in. : 

Forward compartment 

Center compartment 

Rearward compartment 

Midsection 

Maximum external diameter of tank, in 

Xg a , percent chord . 

x cg , percent chord: 

All tanks empty 

Forward tank full 

Rearward tank full 

All tanks full 

kfc, lb/ft 

kg,, ft -lb /rad 1 ar 

m 1 , lb-sec 2 /ft: 

All tanks empty 

All tanks full 

m, lb-sec^/ft: 

All tanks empty 

All tanks full 

T~ , ft-lb-sec 2 : 


0.5 

2 

65-010 


All tanks empty . . . 
All tanks full (solid) 
f h , cps: 

All tanks empty . . . 
All tanks full . . . 
f a , cps: 

All tanks empty . . . 
All tanks full . . . 

g h,0 

8a, 0 * * * 


2.205 

3.33 

4.55 


43-5 

36.4 

52.0 

42.8 

1,955 





0.294 

0.497 

0.1203 

0.323 

0.00862 

0.0209 

15.0 

9-97 


15.8 

10.35 

0.013 
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TABLE II.- GEOMETRIC PROPERTIES OF EXTERNAL TANK 

OF MODEL B 


(a) Tank ordinates 


Distance from 
nose of tank, 
in. 

x, 

ft 

0 

0 

1.09 

.0658 

2.03 

.1016 

' 3-41 

.1358 

4.54 

.1542 

5.88 

.1684 

8.58 

• 1750 

9-875 

-1770 

18.33 

.1770 

20.47 

.1716 

23.77 

-1275 

27.18 

.0642 

30.3 

0 


(b) Volumetric properties 











TABIE III.- MASS MCMEHE-OP-IMERTIA. BATA FOR MODEL A 


Distribution 


Ratios of mass moments 
of inertia 


Emptying Fullness, 

sequence percent 


m 


liffil 






kp . a. 

0 

Effective 

ES9 

2.300 

2.423 

2.232 

2.352 

2.155 

2.280 

2.155 

2.231 

2.106 

2.181 

a 1 . 6 l 4 

s 

1.591 

1.000 

1.000 

| 2.300 

2.375 

1 2.271 

2.324 

1 1 * 759 

1.733 

1.140 

1.143 

1.072 

1.071 

"•s * 5 

1.832 

1.217 

1.242 

1.159 

1.170 

1.092 

1.099 

1.063 

1.049 

1.755 

I.690 

2.203 

2.252 

1.188 

1.192 






























































































£ 

xabu it.- Btnanm. Honan 

{a) Hotel A 
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IABIE IV.- EXPHHMEETEAL RESUIZPB - Continued 
(b) Model B - Fuel-loaded configuration 


= 2,015 lb/ft j Sq. = 543.5 ft-lb/radlan] 


Will T npnn r 

Frequencies 

Flutter 

amplitudes 

Damping 

coefficients 

Ib-sec^'W*' 

v f , 

fps 

D 

percent 

cps 

cps 

cps 

translation; 

in. 

Pitch, 

deg 

Sh 

Ba 

0 

11.2 

16.7 

12.9 

1 

1 

1 

— 

0.009 

0.018 

0.002448 

jjgp 

5.36 




U.4 

0.19 

I 

0.029 

0.072 


178 

4.97 




11.4 

.20 

STM 

.0295 

• 073 

.002461 

181 

5.06 

25 

ID. 35 

13-7 

11-3 

.23 

mm 

.0305 

.077 

.002457 

188 

5-30 




11.4 

.29 

m 

.034 

.080 

.002458 

191 

5-335 




10.2 

0.24 

2.0 

0.027 

0.078 

0.002465 

152 

4.74 

50 

9.64 

11.4 

10.2 

• 30 

2.1 

.030 

.079 

.002462 

161 

5-025 




10.2 

-33 

2.3 

.0345 

.082 

.002465 

168 

5-24 




9.21 

0.27 

2.0 

0.027 

0.074 

0.002475 

132 

4.56 

75 

9.06 

10.1 

9.21 

•35 

2.1 

.032 

.0745 

.002472 

l4l 

4.88 




9.21 

-38 

2.2 

.039 

•0755 

.002471 

l48 

5-12 




8.91 

0.155 

0.8 

0.026 

0.060 

0.002525 

120 

4.29 




8.77 

•235 

1.1 

.026 

.0615 

.002526 

129 

4.68 

ftc: 

8 70 

a.li^ 

a. 77 

■ 255 

1.0 

.026 

.061 

.002520 

133 

4.83 



7* v 

8.78 

• 27 

.95 

.026 

.061 

.002519 

143 

5.18 




8.72 

•32 

.8 

.029 

.060 

.002513 

15? 

5-805 




8.65 

.29 

.65 

.027 

.060 

.002509 

174 

6.405 




8.82 

0.10 

0-55 

0.026 

0.050 

0.002519 

119 

4.295 




8.78 

.12 

.6 

.026 

.050 

.002518 

125 

4.53 

on 

ft, ^2 

Q.OQ 

8.64 

.145 

.45 

.026 

.050 

.002514 

143 

5.27 




8.62 

■ 15 

.4 

.026 

.0495 

.002509 

160 

5-91 




8.55 

■ 15 

• 3 

.026 

-.049 

.002503 

176 

6.56 




8.57 

.14 

.2 

.026 

.0485 

.002494 

202 

7.505 




8.54 

0.06 

0.25 

0.026 

0.038 

0.002519 

101 

3.765 




8.57 

.09 

• 3 

.026 

.038 

.002514 

129 

4-79 

95 

8.45 

8.26 

8.55 

• 06 

■15 

.026 

.038 

.002510 

143 

5-325 




c 8.4o 

.02 

0 

.026 

• 0375 

.002506 

164 

6.22 




°8.35 

.02 

0 

.026. 

• 0375 

.002497 

191 

7-30 

100 

8.36 

7.22 


— 

— 

0.027 

0.026 

0.002393 

d 324 

Q 


^Reference flutter speed V f 0 (see figs. 9 and 11). 

c Bendlng freq uen cies, flutter not discernible from oscillograph, records. 
^Maximum tunnel airspeed, no flutter. 






















































































TABIE IV— EXPERIMENTAL RESULTS - Concluded 


& 


(c) Model B - Effect of tank center of gravity on tank-empty configuration 
jkjj ■ 3,540 lb /ft; k^ a 204.1 ft-lb/radlanj 


TVink 

center-of- 
gravlty 
position, 
percent wing 
chord 

X 

a 

Frequencies 

Damping 

coefficients 

P, 

lb-sec 2 /!!^* - 

V 

fps 

1 

k 

f h> 

cps 

cps 

cps 



50 

0.0235 

14.2 

12.4 


0.008 

0.020 

0.002378 

e 275 

MMMM 

40 

.0767 

14.2 

12.6 

13.7 

.009 

.022 

.002284 

295 

6.87 

50 

.128 

14.1 

12.4 

15.3 

.015 

.021 

.002315 

233 

5.56 

60 

.182 

14.1 

12.15 

13.0 

.008 

.022 

.002404 

256 

5.78 


divergence (based on q. =* 90 lb/ft 2 and standard conditions), no flutter. 
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IABIE 7.- ANALYTICAL RE SUITS 
(a) Model A 


Emptying Fullness, 
sequence percent 




^ ; 


I' • " "V-'v^N'r 






- • I 


Without damping 

Yr* 

fps 

cps 

i 

125 

10.8 

3.63 

155 

10.9 

3.89 

142 

11.0 

4.n 

148 

11.6 

4.04 

144 

12.3 

3.73 

127 

15.3 

3.05 

125 

14.5 

2.775 

127 

n.o 

5.66 

126 

11.6 

3.46 

104 

12.5 

2.66 

109 

13.1 

2.64 

11B 

15.7 

2.74 

no 

n .2 

3.15 

122.5 

11.7 

3.33 

127 

12.0 

3.36 

139 

12.7 

3.47 

125 

15.6 

2.94 

126 

11.7 

5.42 

136 

11.9 

5.64 

115.5 

12.3 

2.98 


With damping 
(& = Sb . °r go ) 


0.0255 


• 0195 


.022 


0.051 


0.024 


0.021 


.059 


.022 


.050 


0.0515 


.024 


0.055 


0.054 


0.022 


0.029 


V 


1 

fps 

cps 

k 

127 

10.6 

3.83 

136 

10.7 

4.05 

145 

10.8 

4.27 

149 

11.5 

4.38 

147 

11-9 

3.95 

l4o 

12.7 

3.50 

128 

14.2 

2.89 

131 

10.8 

3.85 

131 

11.3 

3.69 

118 

12.0 

3-13 

115 

13.1 

2.75 

123 

13.55 

2.89 

117 

10.9 

3.40 

127 

H.7 

3.44 

131 

12.0 

3.47 

l4l 

12.65 

3.55 

130 

13.4 

3.10 

132 

11.4 

3.67 

139 

11.7 

3.77 

122 

12.3 

3.14 
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TABLE V.- ANALYTICAL RESULTS - Concluded. 


(b) Model B - Effect of Including aerodynamic forces 
on external tank containing fuel 


Fullness, 

percent 

Without tank 
aerodynamic 
forces 

With tank 
aerodynamic 
forces 

V f> 

fps 

f fJ 

cps 

1 

k 

v f> 

fps 

ff> 

cps 

1 

k 

0 

223 

13.0 

5.46 

209 

13.1 

5.05 

25' 

198 

ll. 4-5 

5.50 

182 

11.4 

5.07 

50 

163 

10.5 

4.94 

152 

10.6 

4.57 

75 

129.5 

9.63 

4.28 

121.5 

9.60 

4.03 

85 

120 

8.80 

4.3^ 

112 

9.07 

3.92 

90 

95 

100 

104 

8.85 

3.74 

94 

8. 87 

3.36 


(c) Model B - Effect of tank position on flutter and 
divergence of tank-empty configuration 



Tank 

center-of- 

gravity 

position, 

percent 

wing 

chord 


^a,0 

Without tank 
aerodynamic . . 
forces 

With tank 
aerodynamic 
forces 

fps 

Vf, 

fps 

*f* 

cps 

1 

k 

•\ CD 

•s CD 

*fi 

cps 

1 

k 

30 

0.0233 

1.019 

292 




269 




40 

.0767 

1.019 

305 

ON 

'iS 

13.6 

6.30 

278 

244 

13.5 

5.75 

50 

.128 

1.040 

300 

196 

13.4 

4.65 

274.5 

198 

13.5 

4.67 

60 

.182 

H 

• 

0 

• — 3 

0 

294 

183 

13.5 

4.32 

268 

202 

13.4 

4.80 












TANK 

Figure 2- - Schematic ylev of flutter model B and apparatus (fig. 2 of ref. 
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Figure 10.- Effect of fluid on flutter amplitudes of model B without 
haffles. The dashed curves indicate extrapolations. Open symbols 
denote initial flutter points and solid symbols denote limited- 
amplitude flutter. 




Figure 11.- Comparison of experimental and calculated flutter and 
divergence speeds for model B without fuel. 
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(b) 25 percent full; 

^ = 0.75. 

<% 


(c) 50 percent full; 

^ = o.e*. 


Figure 12.- Correlation of theoretical flutter solutions with 
experiment; including damping and aerodynamic tank forces, 
for model B. 
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Figure 12.- Concluded. 
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